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Abstract: We develop a full theoretical analysis of the nonlinear inter- 
actions of the two polarizations of a waveguide by means of a vectorial 
model of pulse propagation which applies to high index subwavelength 
^si . waveguides. In such waveguides there is an anisotropy in the nonlinear 

behavior of the two polarizations that originates entirely from the waveguide 
"^ . structure, and leads to switching properties. We determine the stability prop- 

O ' erties of the steady state solutions by means of a Lagrangian formulation. 

We find all static solutions of the nonlinear system, including those that 
^H' are periodic with respect to the optical fiber length as well as nonperiodic 

soliton solutions, and analyze these solutions by means of a Hamiltonian 



o 



Y^ ' formulation. We discuss in particular the switching solutions which lie near 

the unstable steady states, since they lead to self -polarization flipping which 
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can in principle be employed to construct fast optical switches and optical 



H . logic gates. 



©2013 Optical Society of America 

OCIS codes: (190.4360) Nonlinear optics, devices; (190.4370) Nonlinear optics, fibers; 
(190.3270) Kerr effect; (130.4310) Integrated optics (Nonlinear); (060.4005) Microstructured 
fibers; (060.5530) Pulse propagation and temporal solitons 



\^ \ References and links 

1 . R. H. Stolen, J. Botineau, and A. Ashkin, "Intensity discrimination of optical pulses with birefringent fibers," 
Opt. Lett. 7, 512-514 (1982). 

2. F. Matera and S. Wabnitz, "Nonlinear polarization evolution and instability in a twisted birefringent fiber," Opt. 
Lett. 11, 467^69(1986). 

3. H. G. Winful, "Polarization instabilities in birefringent nonlinear media: application to fiber-optic devices," Opt. 
Lett. 11, 33-35 (1986). 

4. C. Menyuk, "Nonlinear pulse propagation in birefringent optical fibers," Quantum Electronics, IEEE Journal of 
K> ■ 23, 174 - 176 (1987). 

5_J ^ 5. S. F. Feldman. D. A. Weinberger, and H. G. Winful, "Polarization instability in a twisted birefringent optical 

j^ [ fiber," J. Opt. Soc. Am. B 10, 1 191- 1201 (1993). 

" " " 6. G. Millot, E. Seve, and S. Wabnitz, "Polarization symmetry breaking and pulse train generation from the modu- 

lation of light waves," Phys. Rev. Lett. 79. 661-664 (1997). 

7. G. Millot, E. Seve, S. Wabnitz, and M. Haeltennan, "Dark-soliton-like pulse-train generation from induced mod- 
ulational polarization instability in a birefringent fiber," Opt. Lett. 23, 511-513 (1998). 

8. S. Pitois, G. Millot, and S. Wabnitz, "Polarization domain wall solitons with counterpropagating laser beams," 
Phys. Rev. Lett. 81, 1409-1412 (1998). 

9. S. Pitois, G. Millot, and S. Wabnitz, "Nonlinear polarization dynamics of counteipropagating waves in an 
isotropic optical fiber: theory and experiments," J. Opt. Soc. Am. B 18, 432^43 (2001). 

10. S. Wabnitz, "Polarization domain wall solitons in elliptically birefringent optical fibers," PIERS ONLINE 5, 
621-624 (2009). 



11. V. V. Kozlov and S. Wabnitz, "Theoretical study of polarization attraction in high-birefringence and spun fibers," 
Opt. Lett. 35, 3949-3951 (2010). 

12. J. Fatome, S. Pitois, P. Morin, and G. Millot, "Observation of light-by-light polaiization control and stabilization 
in optical fibre for telecommunication applications," Opt. Express 18, 15311-15317 (2010). 

13. V. V. Kozlov, J. Nuiio, and S. Wabnitz, "Theory of lossless polarization attraction in telecommunication fibers," 
J. Opt. Soc. Am. B 28, 100-108 (201 1). 

14. V. E. Zakharov and A. V. Mikhailov, "Polarization domains in nonlinear optics," JETP Lett. 45, 349-352 (1987). 

15. S. Pitois, A. Picozzi, G. Millot, H. R. Jauslin, and M. Haelterman, "Polarization and modal attractors in conser- 
vative counterpropagating four-wave interaction," Europhysics Letters 70, 88-94 (2005). 

16. S. Pitois, J. Fatome, and G. Millot, "Polaiization attraction using counter-propagating waves in optical fiber at 
telecommunication wavelengths," Opt. Express 16, 6646-6651 (2008). 

17. S. Wabnitz, "Cross-polarization modulation domain wall solitons for WDM signals in birefringent optical fibers," 
Photonics Technology Letters, IEEE 21, 875 -877 (2009). 

18. E. Seve, G. Millot, S. Trillo, and S. Wabnitz, "Large-signal enhanced frequency conversion in birefringent optical 
fibers: theory and experiments," J. Opt. Soc. Am. B 15, 2537-2551 (1998). 

19. G. Gregori and S. Wabnitz, "New exact solutions and bifurcations in the spatial distribution of polarization in 
third-order nonlinear optical interactions," Phys. Rev. Lett. 56, 600-603 (1986). 

20. S. M. Jensen, "The nonhnear coherent coupler," IEEE Journal of Quantum Electronics 18, 1580-1583 (1982). 

21. C. M. de Sterke and J. E. Sipe, "Polarization instability in a waveguide geometry," Opt. Lett. 16, 202-204 (1991). 

22. Y. Wang and W. Wang, "Nonlinear optical pulse coupling dynamics," J. Lightwave Technol. 24, 2458 (2006). 

23. Y. S. Kivshar, "Switching dynamics of solitons in fiber directional couplers," Opt. Lett. 18, 7-9 (1993). 

24. D. C. Hutchings, J. S. Aitchison, and J. M. Amold, "Nonhnear refractive coupling and vector solitons in 
anisotropic cubic media," J. Opt. Soc. Am. B 14, 869-879 (1997). 

25. H. S. Eisenberg, Y Silberberg, R. Morandotti, A. R. Boyd, and J. S. Aitchison, "Discrete spatial optical solitons 
in waveguide arrays," Phys. Rev. Lett. 81, 3383-3386 (1998). 

26. U. Peschel, R. Morandotti, J. M. Amold, J. S. Aitchison, H. S. Eisenberg, Y Silberberg, T. Pertsch, and F. Lederer, 
"Optical discrete solitons in waveguide arrays. 2. Dynamic properties," J. Opt. Soc. Am. B 19, 2637-2644 (2002). 

27. K. R. Khan, T. X. Wu, D. N. Christodoulides, and G. I. Stegeman, "Soliton switching and multi-frequency 
generation in a nonlinear photonic crystal fiber coupler," Opt. Express 16, 9417-9428 (2008). 

28. C. C. Yang, "All-optical ultrafast logic gates that use asymmetric nonhnear directional couplers," Opt. Lett. 16, 
1641-1643 (1991). 

29. T. Fujisawa and M. Koshiba, "All-optical logic gates based on nonlinear slot-waveguide couplers," J. Opt. Soc. 
Am. B 23, 684-691 (2006). 

30. W. Fraga, J. Menezes, M. da Silva, C. Sobrinho, and A. Sombra, "All optical logic gates based on an asymmetric 
nonlinear directional coupler," Optics Communications 262, 32 - 37 (2006). 

31. D. C. Hutchings and B. S. Wherrett, "Theory of the anisotropy of ultrafast nonlinear refraction in zinc-blende 
semiconductors," Phys. Rev. B 52, 8150-8159 (1995). 

32. G. P. Agrawal, Nonlinear Fiber Optics (Academic Press, 2007). 

33. S. Afshar V. and T. M. Monro, "A full vectorial model for pulse propagation in emerging waveguides with 
subwavelength structures part I: KeiT nonhnearity," Opt. Express 17, 2298-2318 (2009). 

34. J. B. Driscoll, X. Liu, S. Yasseri, 1. Hsieh, J. I. Dadap, and R. M. Osgood, "Large longitudinal electric fields (E-) 
in silicon nanowire waveguides," Opt. Express 17, 2797-2804 (2009). 

35. B. A. Daniel and G. P. Agrawal, "Vectorial nonlinear propagation in silicon nanowire waveguides: polarization 
effects," J. Opt. Soc. Am. B 27, 956-965 (2010). 

36. V. R. Almeida, Q. Xu, C. A. Barrios, and M. Lipson, "Guiding and confining light in void nanostracture," Optics 
Letters 29, 1209(2004). 

37. O. Boyraz, P. Koonath, V. Raghunathan, and B. Jalali, "All optical switching and continuum generation in silicon 
waveguides," Opt. Express 12, 4094 (2004). 

38. C. Koos, P. Von'eau, T. Vallaitis, P. Dumon, W. Bogaerts, R. Baets, B. Esembeson, I. Biaggio, T Michinobu, 
F. Diederich, W. Freude, and J. Leuthold, "All-optical high-speed signal processing with silicon-organic hybrid 
slot waveguides," Nature Photonics 3, 216 (2009). 

39. W. Astar, J. B. Driscoll, X. Liu, J. 1. Dadap, W. M. J. Green, Y. A. Vlasov, G. M. Carter, and R. M. Osgood, 
"Tunable wavelength conversion by XPM in a silicon nanowire, and the potential for XPM-multicasting," J. 
Lightwave Technol. 28, 2499-251 1 (2010). 

40. R. K. W. Lau, M. Menard, Y. Okawachi, M. A. Foster, A. C. Tumer-Foster, R. Salem, M. Lipson, and A. L. 
Gaeta, "Continuous-wave mid-infrared frequency conversion in silicon nanowaveguides," Opt. Lett. 36, 1263- 
1265(2011). 

41. M. Pelusi, F Luan, T. D. Vo, M. R. E. Lamont, S. J. Madden, D. A. Bulla, D.-Y Choi, B. Luther-Davis, and B. J. 
Eggleton, "Photonic-chip-based radio-frequency spectrum analyser with terahertz bandwidth," Nature Photonics 
3, 139(2009). 

42. X. Gai, T. Han, A. Prasad, S. Madden, D.-Y. Choi, R. Wang, D. Bulla, and B. Luther-Davies, "Progress in optical 
waveguides fabricated from chalcogenide glasses," Opt. Express 18, 26635-26646 (2010). 



43. B. J. Eggleton, B. Luther-Davies and K. Richardson, "Chalcogenide photonics," Nature Photonics 5, 141-148 
(2011). 

44. P. Petropoulos, T. M. Monro, W. Belardi, K. Fumsawa, J. H. Lee, and D. J. Richardson, "2R-regenerative all- 
optical switch based on a highly nonlinear holey fiber." Opt. Lett. 26, 1233-1235 (2001). 

45. H. Ebendorff-Heidepriem, P. Petropoulos, S. Asimakis, V. Finazzi, R. C. Moore, K. Frampton, F. Koizumi, D. J. 
Richardson, and T. M. Monro, "Bismuth glass holey fibers with high nonlinearity," Opt. Express 12, 5082 (2004). 

46. S. Afshar V., W. Q. Zhang, H. Ebendorff-Heidepriem, and T. M. Monro, "Small core optical waveguides are 
more nonhnear than expected: experimental confirmation," Opt. Lett. 34, 3577-3579 (2009). 

47. G. Qin, X. Yan, C. Kito, M. Liao, T. Suzuki, A. Mori, and Y. Ohishi, "Highly nonlinear tellurite microstruc- 
tured fibers for broadband wavelength conversion and flattened supercontinuum generation," Joumal of Applied 
Physics 107, 043108(2010). 

48. F. Poletti, X. Feng, G. M. Ponzo, M. N. Petrovich, W. H. Loh, and D. J. Richardson, "All-solid highly nonlinear 
singlemode fibers with a tailored dispersion profile," Opt. Express 19, 66-80 (201 1). 

49. M. D. Turner, T. M. Monro, and S. Afshar V., "A full vectorial model for pulse propagation in emerging 
waveguides with subwavelength structures part IF Stimulated Raman scattering," Opt. Express 17, 1 1565-1 1581 
(2009). 

50. W. Q. Zhang, M. A. Lohe, T. M. Monro, and S. Afshar V., "Nonlinear polarization bistability in optical 
nanowires," Opt. Lett. 36, 588-590 (2011). 

51. S. Afshar V., W. Q. Zhang, and T. M. Monro, "Structurally-based nonlinear birefringence in waveguides with 
subwavelength stnictures and high index materials," in "ACOFT 2009 Proceeding," (Australian Optical Society, 
2009), pp. 374-375. 

52. E. Magi, L. Fu, H. Nguyen, M. Lamont, D. Yeom. and B. Eggleton, " Enhanced Kerr nonhnearity in sub- 
wavelength diameter As2Se3 chalcogenide fiber tapers," Opt. Express 15, 10324-10329 (2007). 

53. S. Coleman, "Classical lumps and their quantum descendents", in New phenomena in subnuclear physics, Ed. A. 
Zichichi (New York, 1977), pp. 185-264. 

54. Y S. Kivshar and G. P. Agrawal, "Optical soHtons: from fibers to photonic crystals" (Academic Press, 2003) 

55. W. Q. Zhang, M. A. Lohe, T. M. Monro, and S. Afshar V., "Nonlinear self-flipping of polaiization states in 
asymmetric waveguides", to be published. 

56. W. Q. Zhang, M. A. Lohe, T. M. Monro, and S. Afshar V., "Nonlinear polarization self-flipping and optical 
switching," in Proceedings of the International Quantum Electronics Conference and Conference on Lasers and 
Electro-Optics Pacific Rim 2011, (Optical Society of America, 2011), paper C370. 

57. I. S. Gradshteyn and I. M. Ryzhik, "Table of Integrals, Series and Products" (Academic Press, 1965). 

58. W. Q. Zhang, M. A. Lohe, T M. Monro, and S. Afshar V., " New regimes of polarization bistability in linear 
birefringent waveguides and optical logic gates," in Nonlinear Photonics, OSA Technical Digest (CD) (Optical 
Society of America, 2010), paper NThD4. 



1. Introduction 

The Kerr nonlinear interaction of the two polarizations of the propagating modes of a wave- 
guide leads to a host of physical effects that are significant from both fundamental and applica- 
tion points of view. Here, we develop a model of nonlinear interactions of the two polarizations 
using full vectorial nonlinear pulse propagation equations, with which we analyze the nonlinear 
interactions in the emerging class of subwavelength and high index optical waveguides. Based 
on this model we predict an anisotropy that originates solely from the waveguide structure, and 
which leads to switching states that can in principal be used to construct optical devices such 
as switches or logical gates. We derive the underlying nonlinear Schrodinger equations of the 
vectorial model with explicit integral expressions for the nonlinear coefficients. We analyze 
solutions of these nonlinear pulse propagation equations and the associated switching states 
by means of a Lagrangian formulation, which enables us to determine stability properties of 
the steady states; this formulation provides a global view of all solutions and their properties 
by means of the potential function and leads, for example, to the emergence of kink solitons 
as solutions to the model equations. We also use a Hamiltonian formalism in order to identify 
periodic and soUtonic trajectories, including solutions that allow polarization flipping, and find 
conditions under which the unstable states and associated switching solutions are experimen- 
tally accessible. 

The nonlinear interactions of the two polarizations of the propagating modes of a waveguide 
have been studied extensively over the last 30 years ifTI- lfTsl . Different aspects of the inter- 



actions have been investigated, for example Stolen et al. |r| used the induced nonlinear phase 
difference between the two polarizations to discriminate between high and low power pulses. In 
the context of counterpropagating waves, the nonlinear interactions have been shown to lead to 
polarization domain wall solitons, lISl- lfTOl . lfT4l which are described as kink solitons represent- 
ing a polarization switching between different domains with orthogonal polarization states. The 
nonlinear interactions can also lead to polarization attraction ||9], lfTTl - llT3l . llT5l[T6l where the 
state of the polarization of a signal is attracted towards that of a pump beam. For twisted bire- 
fringent optical fibers, polarization instabiUty (|2] |5] and polarization domain wall solitons ifTTl 
have been reported. The nonlinear interactions also induce modulation instability which results 
in dark-soliton-like pulse-train generation [:6. 7|. Large-signal enhanced frequency conversion 
ifTSl . cross-polarization modulation for WDM signals ifTOl . and polarization instability IS) have 
also been reported and attributed to nonlinear polarization interactions. Stability behavior has 
been studied in anisotropic crystals lfT9l . 

The nonlinear interactions of the two polarizations can also be studied in the context of either 
nonlinear coherent coupling or nonlinear directional coupling in which the amplitudes of two 
or more electric fields, either the two polarizations of a propagating mode of a waveguide or 
different modes of different waveguides, couple to each other through linear and nonlinear 
effects II20I - II22II . Nonlinear directional coupling is relevant to ultrafast all-optical switching, 
such as soliton switching ll23]| - ll27l and all-optical logic gates ll28l - ll30l . The interaction of 
ultrafast beams, with different frequencies and polarizations, in anisotropic media has also been 
studied and the conditions for polarization stability have been identified ll24ll3Tl . 

In previous work ( ll32l . Chapter 6), the nonlinear interactions of the two polarizations are 
described by two coupled Schrodinger equations. These equations employ the weak guidance 
approximation, which assumes that the propagating modes of the two polarizations of the wave- 
guide are purely transverse and orthogonal to each other within the transverse x^y plane, per- 
pendicular to the direction of propagation z- Based on this, the electric fields are written as 

E,-(x,y,z,0=A,(z,f)e;(jc,3;), /-=1,2, (1) 

where A,(z,f) are the amplitudes of the two polarizations, with t\{x,y) . e2(-*:,}') = 
e\{x,y)e2{x,y)^ .y = 0, where e\{x,y),e2{x,y) are the transverse distributions of the two po- 
larizations, x,y are unit vectors along the x and y directions, and it is understood that fast 
oscillatory terms of the form exp(— /OJf ± j3,z) are to be included for the polarization fields. 
The weak guidance approximation also assumes that the Kerr nonUnear coefficients for the self 
phase modulation of the two polarizations are equal because their corresponding mode effective 
areas are equal ll32l . We refer here to models of nonlinear pulse propagation based on the weak 
guidance approximation simply as "scalar" models, since these models consider only purely 
transverse modes for the two polarizations. 

The weak guidance approximation works well only for waveguides with low index con- 
trast materials, and large dimension structure compared to the operating wavelength. This ap- 
proximation is, however, no longer appropriate for high index contrast subwavelength scale 
waveguides (HIS-WGs) lf33l - ll35l . These waveguides have recently attracted significant inter- 
est mainly due to their extreme nonlinearity and possible applications for all optical photonic- 
chip devices. Examples include silicon, chalcogenide, or soft glass optical waveguides, which 
have formed the base for three active field of studies: silicon photonics ll36l - P0l . chalcogenide 
photonics BTl - ll43l . and soft glass microstructured photonic devices B4l - ll48l . 

In order to address the limitations of the scalar models in describing nonlinear processes 
in HIS-WGs, we have developed in ||33l a full vectorial nonlinear pulse propagation model. 
Important features of this model are: (1) the propagating modes of the waveguide are not, in 
general, transverse and have large z components and, (2) the orthogonality condition of different 



polarizations over the cross section of the waveguide is given by /ci {x,y) x h2{x,y) . z dA = 
0, rather than simply ei(x,3') .e2{x,y) = as in the scalar models. These aspects lead to an 
improved understanding of many nonlinear effects in HIS-WGs; it was predicted in ||33l . for 
example, that within the vectorial model the Kerr effective nonlinear coefficients of HIS-WGs 
have higher values than those predicted by the scalar models due to the contribution of the z- 
component of the electric field, as later confirmed experimentally ||46l. Similarly, it was also 
predicted that modal Raman gain of HIS-WGs should be higher than expected from the scalar 
model ED. 

Here, we extend the vectorial model to investigate the nonlinear interaction of the two po- 
larizations of a guided mode. The full vectorial model leads to an induced anisotropy on the 
dynamics of the nonlinear interaction of the two polarizations |50|, which we refer to as struc- 
turally induced anisotropy, in order to differentiate this anisotropy from others, such as those 
for which the anisotropy originates from isotropic materials. The origin of the anisotropy is the 
structure of the waveguide rather than the waveguide material. 

The origin of this anisotropy in subwavelength and high index contrast waveguides has also 
been reported by Daniel and Agrawal [35], who considered nonlinear interactions of the two 
polarizations in a silicon rectangular nanowire including the effect of free carriers. In their 
analysis, however, they ignore the coherent coupling of the two polarizations, considering the 
dynamics of the Stokes parameters only for a specific waveguide and ignore the linear phase. 

This anisotropy in turn leads to a new parameter space in which the interaction of the two 
polarizations shows switching behavior, which is a feature of the vectorial model not accessible 
through the scalar model with the underlying weak guidance approximation. We also show that 
the resulting system of nonlinear equations, for the static case, can be solved analytically. Due to 
the underlying similarity of the nonlinear interaction of the two polarizations and the nonlinear 
directional coupling of two waveguides, the anisotropy discussed here can be also applied to 
the case of nonlinear directional coupling, in which the two waveguides have different effective 
nonlinear coefficients for the propagating modes. 

This work develops and expands on results reported for the first time in ll50ll5Tl . in particular 
we derive here (in Section |2) the equations that describe the nonlinear interactions of the two 
polarizations within the framework of the vectorial model, with explicit integral expressions for 
the nonlinear coefficients. In Section[3]we determine properties of the static solutions, classify 
the steady state solutions, and determine their stability using a Lagrangian formalism. We also 
discuss a Hamiltonian approach and how the phase space portrait provides a complete picture 
of the trajectories of the system, including the periodic and solitonic solutions (Section lSTSl l. We 
derive analytical periodic solutions by direct integration of the system of equations in SectionlH 
and then discuss switching solutions and their properties. We relegate to the Appendix a math- 
ematical analysis of the exact soliton solutions, which are relevant to the switching solutions, 
with concluding remarks in Section |5] 

2. Nonlinear differential equations of the model 

In the vectorial model the nonlinear pulse propagation of different modes of a waveguide is 
described by the equations: 

dAy ~ i"-^pl"^ d"Ay 
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where /i,V = 1,2 with jj. ^v, and Ai{z,t),A2{z,t) are the amplitudes of the two orthogonal 
polarizations. These equations follow from the analysis in ll33l . by combining Eqs. (23,32) 
of II33I . but without the shock term. The linear birefringence is defined by AjSy^ — — Aj3^v = 
j3v — Pn and the 7 coefficients are given by 
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Here we use the notation (ey)^ = ey .Cy, |eyp = Cy .Cy and |ey| = (cy .ey)(ey .Cy), together 

I |2 

with |ey .Cul = (ey.e?()(ey .e^). In these equations ei(jt:,y), 62 (a:,}') are the modal fields of the 
two orthogonal polarizations, k = In/X is the propagation constant in vacuum, and yy, y^y, 

Yuv I 7ai V ' Tfiv 1 Yfiv ^rc the effective nonlinear coefficients representing, respectively, self phase 
modulation, cross phase modulation, and coherent coupling of the two polarizations, and 
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is the normalization parameter 

The coupled equations (|2|i describe the full vectorial nonhnear interaction of the two po- 
larizations. There are two fundamental differences between these equations and the typical 
scalar coupled Schrodinger equations (see for example Chapter 6 in 11321 ). Firstly, the addi- 
tional terms A*Ay,Afi \Av\ ,A^ IA^I on the right hand side of Eq. dU represent interactions 
between the two polarizations. These do not appear in the scalar model since the effective non- 
linear coefficients associated with these terms, 7|iy , 7^v, 7^v ^s given in Eqs. ( 161718b . contain 
factors such as e^ .Cy which are zero in the scalar model, since the modes are assumed to be 
purely transverse. All possible third power combinations of the two polarization fields, namely 
|Ay|^Av,|A^| Av,A'iAl,,A*Al,An\Ay\^ andA^, \Ah\ occur on the right hand side of Eq. ©, 
due to the z-component of the modal fields. Secondly, in all effective nonlinear coefficients 
given by Eqs. (l3][8]l, the modal fields e and h have both transverse and longitudinal compo- 
nents, unlike the scalar model in which modal fields have only transverse components. The 
terms containing nonzero e^ . Cy provide a mechanism for the interaction of the two polariza- 
tions since they allow for exchange of power between the two modes through the z-components 
of their fields. The last term on the right hand side of Eq. (|2), for example, indicates a coupUng 
of power into a polarization, even if initially no power is coupled into that polarization. 

Although the terms on the right hand side of Eq. (|2]i that contain e^j . Cy are nonzero, they 
are generally significantly smaller than the remaining terms and are therefore neglected in the 



following; further investigation of the effects of these terms, and a discussion of their physical 
significance, will be presented elsewhere. The focus of this paper is to investigate the effect of 
the z-components of the fields e and h, which influence the values of the effective coefficients, 
and therefore also the nonlinear interactions of the two polarizations. Hence, from (|2]i, we obtain 
the equations: 

77=1 

These are similar in form to the scalar coupled equations ( ll32l . Section 6.1.2), however, the 
coefficients 7v,7^v,7iiV' given in Eqs. (l3]|5]i, now contain z-components of the electric field, 
through both e and h. In the framework of the scalar model, the weak guidance approximation 
assumes that the effective mode areas of the two polarization modes are equal ||32l , leading to 

71=72 = 37^/2-3/,, (11) 

where we have denoted 7 = 712 = 721, 1^ = 1^2 = ^21- This means that in the scalar model 
there is an isotropy of the nonlinear interaction of the two polarizations; in order to break this 
isotropy, one needs to use either anisotropic waveguide materials or twisted fibers, or else cou- 
ple varying light powers into the two polarizations by using either counter- or co-propagating 
laser beams. The fact that in the vectorial form (IT&i of the coupled equations the 7 values in- 
clude the z-component of the fields, as given by Eqs. (l3]|5]l, means that the equalities ( fTTT ) do 
not hold in general. As an example, see Fig. 1 in (501 which plots 71 , 72, 7c, 1^ for a step-index 
glass-air waveguide with an elliptical cross section; evidently the equalities (fTTT i are not satis- 
fied. One consequence of the vectorial formulation is, as we show in Section [3341 the existence 
of unstable states not present in the scalar formulation. 

3. Static equations 

We find now all solutions of Eq. ( fTOl l for the static case, in which the fields Ai ,A2 are functions 
of z only. We have therefore the two equations 

^ = /(7i|A,|2 + 7|A2p)A,+//A2At.-2'-^'^-^ (12) 

^ = /(72|A2p + 7:|Aip)A2 + //A2A|.2'^/^'-, (13) 

where Aj3 = j5i — ^2- We express the fields Ai ,A2 in polar form according to 

Ai=ViYe'*, A2^^/P2e'^, (14) 

where the powers Pi ,P2 and the phases 0i , 02 are real functions of z- It is convenient to define 
the phase difference A0 and an angle 9 according to 

A0 = (/)!- 02 +zAi3, 0=2A0, (15) 

then upon substitution into Eqs. ( 112113b we obtain the four real equations: 

^ = 2/.PiP2sin0 (16) 

dz 

dPo 

-^ = -2i/.PiP2sin0 (17) 

dz 

dB 

— = 2Aj3+2Pi(7i-7,-i/.cos0)-2P2(72-7c-ycCos0) (18) 

dz 

^ = 7iPi+P2(7+l^cos0). (19) 

dz 



The last equation decouples from the remaining equations, hence we first solve Eqs. (ll6H18b 
for Pi,P2,9 and then determine (pi by integrating ( fT9] l. Eqs. (|16|17t show that Pq = P1+P2 is 
constant in z- We define the dimensionless variables 

Pi Pi 



v = 



Po P1+P2 
and the dimensionless parameters 



2iPoz, (20) 



Aj3 Yc-n , 71+72-27, 
fl = ; ; — , 0= :r-, • (21) 

YcPo Yc ' li. 

In terms of these parameters we obtain the two equations; 

v= — = v(l-v)sin0, (22) 

dx 

dB 

e = — ==-fl + 2^v+(l-2v)cos0. (23) 

Since t takes only positive values, we may regard T as a time variable which is limited in value 
only by the length of the optical fiber and by the value of Pq, and we set the initial values 
I'o ~ v(0), 00 = 0(0) at time T = 0, i.e. at one end of the fiber The general solution depends on 
the initial values vq, Oq and on only two parameters a,b, even though Eqs. ( fT6lfT9] ) depend on 
the five constants -P07 7i , 72, 7 , Yc- 

At the initial time we have Pi , P2 > and so we always choose vq such that < I'o < 1 . It may 
be shown from Eqs. (l22H23b that < v(t) < 1 is then maintained for all T > 0, i.e. the powers 
Pi,P2 remain strictly positive at all later times. The constraint < vq < 1 implies that the initial 
speed 00 is restricted, since it follows from Eq. (l2Jt that 1 1 ^ |a| + 2|fe| + 1 at all times T. 



3.1. Properties of a, b 

Of the two dimensionless parameters a,b, evidently b depends only on the optical fiber param- 
eters, whereas a depends also on the total power Pq, unless Aj3 = 0. For the scalar model, when 
Eqs. (fTTT i are satisfied, we have b — I but generally /? 7^ 1. In this case a set of steady state 
solutions appears (the states JTM discussed in Section [32] below) which for certain values of 
a,b are unstable. For fibers with elliptical cross sections we find that b > 1 and the unstable 
steady states exist provided 1 < a <2b— I. We have not, however, been able to eliminate the 
possibility that b < I for other geometries, and so in the following we also analyze the case 
b < I. The parameter a can be positive or negative depending on the sign of Aj3 and on the 
value Po; when Eqs. ( fTTI ) are satisfied we have a = — 3Aj3/(Po7i ) + 1 and hence a can take large 
positive or negative values for small Pq. 

As an example, we have evaluated b using the definitions Eqs. (|3]|5]l for step-index, air-clad 
glass waveguides with elliptical cross sections where the major/minor axes are denoted x,y. The 
host glass is taken to be chalcogenide with linear and nonlinear refractive indices of n = 2.8 
and«2 = 1.1 X lO^^^m^AV at A = 1.55/im(as in 1521 ). Fig.[ni) shows a contour plot of log 10^ 
as a function of .ic,y. We see, as expected, that b approaches 1 as the waveguide dimensions 
x,y increase towards the operating wavelength. For small core waveguides, however, we find 
b > I with values as large as /? w 200. The parameter a, on the other hand, depends on both the 
structure and the total input power Po. For low input powers, specifically for Pq^. ^ |A/3 1, a can 
take large negative values (for Aj3 > 0) or positive values (for Aj3 < 0) as shown in Fig.llJii)- 
For large values of Pq, however, a approaches the constant C = (72 — yc)lic^ whose contours 
for elliptical core waveguides are shown in Fig. [TJiii); most such waveguides have positive C 
values ranging up to 400, but some, those in the region on the left side of the white curve in 
Fig.fHiii), have negative or small values of C The contour plot for A/3 in Fig. [itiv) shows that 
Aj3 takes a wide range of positive and negative values as jc,y vary. 
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Fig. 1. Contour plots as functions of the elliptical waveguide dimensions x, y of (i) log^^b; 
(ii) a as defined in ^Q for Pq = IW; (iii) C = (}^ - yc)/Yc where C < to the left of the 
white line; (iv) the birefringence Aj3. 



3.2. Steady state solutions 

There are four classes of steady state solutions of Eqs. (I22]23l l, each of which exist only for 
values of a,b within certain limits, as follows: 



cos = 1, V 



provided b y^ 1 and < 



g-l 

2(h-\) 



<l; 



COS0 



-1, V 



fl — 1 



fl + 1 

2(^7+1) 



provided by^ —I and < 



provided \a\ ^ 1; and 



g+l 
2(^+1) 



< 1; 



COS0 = fl, V = 



cosO = —a + 2b, v=l 



(24) 



(25) 



(26) 



(27) 



provided |a — 2fe| ^ 1. 

Of these four classes, ( |26] | and dZTJ i lie on the boundary of the physical region < v < 1, but 
nevertheless influence properties of nearby nontrivial trajectories, and also play a role in soliton 



solutions. The states ( l24l i lie within the physical region only if the parameters {a,b) belong to 
either the red or green region of the a,b plane shown in Fig.|2](i). Similarly the solutions (IZSl l 
satisfy < v < 1 only in the disjoint regions of the a,b plane defined by either 2b+l <a < —I 
or —1 < a <2b+l.U a,blie outside these regions, and also outside the strips given by |a| ^ 1 
and \a — 2b\ ^ 1, there are no steady state solutions. 

For special values of a,b these steady states can coincide, for example if a ~ 1 the solution 
(|26] | coincides with the boundary value of (l24l i. Steady states for values of a,b on the boundary 
of the regions shown in Fig. |2]may need to be considered separately; for example if a = b = 1 
then all steady states are given either by (IZSl l, or else by cos = 1 and any constant v. 

In practice, the values of a,b are determined by the waveguide structure, the propagating 
mode and, in the case of a, the input power Pq, and hence only restricted regions of the a,b 
plane are generally accessible. For example. Fig. [TJi) shows that for the fundamental mode of 
elliptical core fibers we have logjQ^ ^ 0, and so the attainable values of b are limited to b^ 1. 
We nevertheless include the case b <l in our analysis, since this possibility cannot be excluded 
for other fiber geometries. We discuss the accessible regions for the case of unstable steady 
states in Section[331 
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Fig. 2. The a^b plane showing: (i) the regions of existence for the solutions l l24t . either 
1 <a<2b—i (red), or 2fe — 1 < a < 1 (green); (ii) the regions of existence for the unstable 
solutions consisting of M4\ (red), and ( I25t for which 2b + I < a < —1 (orange), together 
with J26l27t for which |a| < 1 or \a -2b\<l (Ught blue). 



3.3. Lagrangian formulation 

We wish to determine the stability properties of each of the four classes of steady state solutions, 
in particular we look for unstable steady states. These are of interest because polarization states 
which lie close to these unstable states are very sensitive to small changes in parameters such as 
the total power Pq, and so can flip abruptly as a function of the optical fiber length z- Although 
we may determine stability properties by investigating perturbations about the constant solu- 
tions, we find it convenient to reformulate the defining equations (I22I23I) as the Euler-Lagrange 
equations of a Lagrangian L which is a function of 0, 0, and depends otherwise only on the 
parameters a,b. This also provides insight into the properties and solutions of these equations, 
and we may then investigate stability by examining the corresponding potential function. From 



(|23] | we have 



+ fl — cos 
2(b-cose) ' 



(28) 



and by substitution into ( l22b we obtain 

2(/7-cos0) 0-sin0 0' + sin0(fl-cos0)(fl-2/7 + cos0) = 0. (29) 

We consider Lagrangians L of the form 

L^T-V = -M{9)0^-V{9) (30) 

where T is the (positive) kinetic energy, V is the potential energy, and the "mass" M is a positive 
function of 9. The equation of motion is 

M{9)9 + ^M'{9)9^ + V'{9)=0, (31) 



and is identical to Eq. ( |29] l provided 

M(e)= ,, ^ ,, , y(0) = -|/,-cos0|--fc^. (32) 

|fc— cos0| |fc— COS0| 

We may therefore investigate all possible solutions 0(t) by analyzing properties of the pe- 
riodic potential V{9); every solution of the system of equations (122123b corresponds to the 
trajectory 0(t) of a particle of variable mass M in the potential V. Steady state solutions are 
zeroes of V'{9), and stability is determined by whether these zeroes are local maximums or 
minimums of V, subject to the constraint that the associated function v should always satisfy 
< V < 1. Trajectories which begin near a local minimum, with a small initial speed 0(0), os- 
cillate periodically with a small amplitude. On the other hand, trajectories which begin near an 
unstable point, i.e. near a local maximum of V, can display periodic oscillations of large ampli- 
tude with abrupt transitions between adjacent local maximums; we refer to these as switching 
solutions (previously bistable solutions i50|) since cosA(p = cos(0/2) switches periodically be- 
tween two distinct values. Soliton trajectories also occur in which the particle moves between 
adjacent local maximums of V, see for example the discussion in 1,53 J , Section 2 and ||54|| for 
properties of solitons in optical fibers. As mentioned in Section [331 soliton trajectories also 
appear as the separatrix in phase plane plots. 

We plot y as a function of and either a or Z? in Fig. [3] showing that V defines a complex 
surface with valleys and peaks which change suddenly as a or Z? are varied. Periodic solutions 
occur for trajectories restricted to a local valley, but there are also unbounded trajectories, in 
which increases or decreases indefinitely, depending on a,b and on whether 0(0) is suffi- 
ciently large. The potential, as a function of and a, has saddle points which indicate that a 
stable solution can become unstable as a is varied; according to the definition (ISTT i we may vary 
a within certain limits by varying the total power Pq. 

For a ~ b the potential is essentially that of the nonlinear pendulum under the influence of 
gravity, namely a simple cosine potential, but with a mass that depends on 0. Provided b > I 
this mass varies between two positive, finite limits. The unstable steady states correspond to 
a pendulum balanced upright, while the switching states (discussed in Section [Hi correspond 
to trajectories which begin with the pendulum positioned near the top, possibly with a small 
initial speed, then swinging rapidly through = 27r to reach the adjacent unstable steady state. 
During this motion cosA0 = cos j flips rapidly between the values ±1. The soliton discussed 
in the Appendix is the trajectory in which the pendulum begins at the unstable upright position 




Fig. 3. The potential V plotted as a function of (i) 6, a for b = 0.8; (ii) d,b for a = 0. 



and, over an infinite time, moves through the stable minimum to the adjoining unstable steady 
state. 

Although both M and V are singular when cos0 = b, which occurs only if \b\ ^ 1, this 
singularity is an artifact of the Lagrangian formulation, as is evident from Eqs. (I22I23I I, which 
have smooth bounded right hand sides for any b. In particular v, which is obtained from Eq. 
( |28] l given 0, is a smooth function of T even if cos 6 = b for some T. 

The energy T + V ^ jM{9) 9^ + V{9) is a constant of the motion. Hence we may integrate 
Eq. (|29ll to obtain 



9'-^{b~ cos 9y + {a-bY + c( 



COS0) 



(33) 



where c is the constant of integration. This constant is determined by first choosing initial values 
vo,9o, where < vq < 1, and then finding 0(0) from Eq. ( |23] | which, from (l33T l evaluated at 
T = 0, fixes c. We may integrate (|33] | to determine 9 as an explicit function of T, expressible in 
terms of elliptic functions, as discussed further in Section [33] 

A Umitation of the Lagrangian formulation is that the constraint < v < 1 is not easily 
implemented. Whereas every solution of the system (I22]23t defines a trajectory 0(t) in the 
Lagrangian system (l29T l, the converse is not true, i.e not all trajectories in this system satisfy 
< V < 1. The initial speed 0(0) must be restricted to only those values allowed by Eq. (l2Jt . in 
which < v(0) < 1, and similarly the constant solutions of Eq. ( l29l l are valid steady states for 
the original equations ( 122123b only in certain regions of the a,b plane. Trajectories which violate 
< I' < 1, while not physical in the context of optical fiber configurations, can nevertheless be 
viewed as acceptable motions of the mechanical system defined by the Lagrangian ( [30] l. We 
investigate an alternative Hamiltonian formulation in terms of v in Section [33] 

3.4. Stability of steady state solutions 

The stability of each of the four classes of steady state solutions in Section [3l2l is determined by 
the sign of V " for that solution; a positive sign implies that the solution lies at a local minimum 
of V and is therefore stable, whereas a negative sign implies that the solution is unstable. 

For the steady states (|24]| we have V" = {a — l){a — 2b+l){b—l)/\b— 1 p and so these states 
are stable for points a,b such that {a — l){a — 2b + l){b — I > 0, shown as the green region in 
the a,b plane in Fig. |2](i), and are unstable in the red region, where 1 < a < 2b— I. For the 
steady states ( |25] | we have V" = (a + 1)(— a + 2b+l){b+ l)/\b+ 1\^ and so these solutions 



are stable for 
in Fig. ID (ii). 



1 < fl < 2^7 + 1 and are unstable in in the orange region 2b+l<a<—l shown 



For the remaining steady states (126127b . for which v = or v = 1, we have Y" = 
— 2sin^ 0/|fl — h\ which in all cases is negative, and so these states are unstable whenever they 
exist. This is consistent with the observation that v(t) cannot attain the values 0, 1 at any time 
T, provided < vq < 1. The regions in the a,h plane where the unstable states exist are shown 
inFig.|2](ii). 

Next, we determine conditions under which the unstable steady state solutions ( |24| ) are ac- 
cessible. For elliptical core step index fibers, for which Z? > 1 as shown in Fig.[TJi), the region 
of instability is indeed accessible and leads to properties such as nonlinear self-polarization 
flipping, discussed in Section ID The region of unstable solutions is given by 1 < a < 2fo— 1, 
equivalently 

Jc + t-l\<^-<li-lc-ic- (34) 

These inequalities specify the possible values, if any, of Pq for which the unstable solutions 
exist for a fixed fiber In order to visualize this region we plot a as a function of Pq in Fig. Ui), 
where a is given by (|2TI) . The boundaries of the unstable region at a = 1 , a = 2/? — 1 are shown 
by the green solid lines. 

First we consider fibers for which 1 < C < 2i> — 1, where C = {yi — Jc^ He takes the value 
shown by the dashed line in Fig.Hfi). Then a has two branches associated with either Aj3 < 
or Aj3 > 0; for the branch corresponding to Aj3 < (the solid blue line), a is large and positive 
for small P{) and asymptotically approaches C for large Po. The intersection of this branch with 
the boundary a=2h —\ determines the minimum power Pmini required in order to access the 
unstable region. In this case, only part of the unstable region corresponding to C < a < 2/? — 1 is 
accessible, as shown by the blue region. For the A/3 > branch (red solid curve) a is large and 
negative for small Pq and asymptotically approaches C for large P^. For this branch, Pq needs 
to be larger than a value Pmin2 • The unstable region is accessible provided 1 < a < C and is a 
subset (red shaded) of the whole unstable solution region. Fig.|4ji) allows one to determine the 
minimum and maximum values of a and the minimum power to access the unstable solution 
region, once A/3 and C are known. For elliptical core fibers these two values are completely 
determined by the dimensions x,y, see Fig. [TJiiijiv) for plots of C and Aj3 . 

Besides fibers for which 1 < C < 2/? — 1, there are the possibilities C >2h —\ orC< 1. 
From Fig. [TliiiJv) one can show that these combinations (with A/3 positive or negative) either 
do not exist, or do not lead to unstable solutions, since the possible values of a do not lie in the 
unstable region 1 < a < 2/? — 1 . In summary, the only elliptical core fibers that allow unstable 
solutions are those with 1 <C <2b~ I with either positive or negative A)3 . The case in which 
Aj3 = is discussed separately in Il55ll56l . 

Based on the above discussion, one can find the minimum power P™'" required to generate 
unstable solutions for elliptical core fibers. Fig. 21ii) plots logio(PQ"'") (where P™" is measured 
in watts) as a function of x,y, where the white region corresponds to fibers for which there are 
no unstable solutions, and the regions below and above the diagonal line correspond to Pmin, 
and Pmin2' respectively, which have been obtained for the two branches of the function a(Po) 
shown inFig.HJi). 

3.5. Hamiltonian function 

Although the Lagrangian formulation in terms of Q is convenient for an analysis of the steady 
states and their stability, and also for a qualitative understanding of all solutions including soli- 
tons, the constraint < v < 1 is more easily implemented by means of a direct formulation 
in terms of v. This automatically eliminates unphysical trajectories for which one of the input 
powers Pi ,P2 is negative. Such a formulation follows by construction of a Hamiltonian function 
which, being conserved, allows us to firstly integrate the nonlinear equations and obtain analyt- 
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Fig. 4. (i) a as a function of Pq for A/3 < (blue solid line) and Aj8 > (red solid line). The 
green lines mark the boundaries of the (red) region of instability in the a, b plane shown in 
Fig-ISi); (ii) contour plot of logio(/o"'") as a function of jc,>', showing the minimum total 
power P™" (in units W) required to access unstable steady states, where they exist. 



ical solutions and, secondly, to interpret physically the possible states of polarizations within 
an optical waveguide from the phase plane contours. Corresponding to the conserved energy 
T +V which follows from (|30] | there is a Hamiltonian function H defined by 



which satisfies 



//(v,0) = -flv + /7v' + v(l-y)cos0 



dH ■ dH 



(35) 



<90' " dv 

Hence as a function of T, H is conserved and takes the constant value Hq ^H{vo, Oq) on any tra- 
jectory. We may investigate all possible solutions, therefore, by analyzing the curves of constant 
Hq in the v, 9 plane. We have 

H() + av — bv^ 



cosO 



v(l-v) ' 



and from Eq. ( |22] | we obtain 



v2 = e(v) 



where Q is the polynomial of 4th degree (provided b^ ^ 1) given by 

Q{v)^v'il-vf-{Ho + av-bvY- 



(36) 

(37) 

(38) 



Since the left hand side of (|37] i is positive, solutions exist only if Q(v) ^0 for v in the interval 
< V < 1. Generally 2(0), 2(1) < but since Q{vo) = v^(l - vo)^ sin^ 0o ^ (as follows from 
Eq. ( |22] |) Q has at least two real zeroes, possibly repeated, and so there is an interval within 
< V < 1 in which Q{v) > 0, and so solutions always exist. If the initial values vq, Oq are such 
that the trajectory begins in a stable steady state, v remains constant for all T > 0, otherwise 
the trajectory is nontrivial. There are two types of nontrivial solutions, periodic and soliton 
solutions. 

We can gain insight into possible solutions by plotting contours of constant H{v,9) in the 
V, 9 plane, which supplies essentially a phase portrait of the system. Solutions for which both 
V, 9 are periodic in T form closed loops, and lie close to a stable steady state, whereas nonpe- 
riodic trajectories lie outside the separatrix which defines soliton solutions, as we discuss in 
the Appendix. Fig. |5] shows two examples in which stable steady states are marked in green. 



and unstable steady states are shown in red or orange. Periodic solutions are evident as closed 
loops surrounding stable steady states, whereas the separatrix marks soliton trajectories which 
connect unstable steady states. Apart from these solitons, all other solutions v,cosO (but not 
necessarily 9) are periodic in T. The switching solutions of particular interest, in which the 
state of polarization inside the waveguide flips between two well-defined states, are those close 
to the separatrix. 





e 

Fig. 5. Contours in the 6,v plane of constant H for (i) fl = l,ii = 4; (ii) a = b = 2, with 
steady states marked by green dots (stable) and red or orange dots (unstable). The separa- 
trix, which identifies the soliton trajectories, is shown in red. 



4. Periodic solutions 

Periodic solutions v of dJTJ l attain both minimum and maximum values, denoted Vmin,Vmax 
respectively, with < Vmin ^ Vmax < 1- Since v = at a maximum or minimum of v, both 
Vmin, Vmax are roots of Q. We can factorize Q as a product of quadratic polynomials. 



Q{v)^^[{b+iy-{a + l)v-HQ\[{b-iy-ia-l)v-Ho\, 



(39) 



and hence explicitly find all roots, and so identify Vmax and Vmij,. We integrate v = \/Q{v) over 
the half-period in which v increases, in order to find T as a function of v, and also the period T: 



du 



VW) 



T-To, 



du 



VW) 



(40) 



where Tq is the time at which v achieves its minimum, i.e. Vmin = v{to). These integrals may 
be evaluated in terms of elliptic integrals of the first kind, see for example the explicit formulas 
in ||57l (Sections 3.145, 3.147). In particular, T is expressible in terms of the complete elliptic 
integral K, and so can be written as an explicit function of a,b,VQ, 9o, i.e. as a function of the 
waveguide parameters and the initial power and phase of the input fields. The precise formulas 
depend on the relative location of the roots of Q. 

Having found v, cos 9 is obtained from Eq. (|36] | and is also periodic in T, as is 9 which is 
obtained from Eq. (l23l l, however 9 itself need not be periodic. Although it is straightforward to 
find V, 9 numerically as functions of T, for specified numerical values of a,b and initial values 



vo, 00. the exact solutions are useful because they display the exact dependence of the solution 
on all parameters, such as the total power Pq; it is not necessary therefore to solve the equations 
numerically for every choice of Pq, rather the exact solution gives the explicit periodic solution 
and the period as known functions of Pq- 

For switching solutions, the phase difference between the two polarization vectors experi- 
ences abrupt phase shifts through n as the light propagates within the waveguide. As a result, 
the state of polarization flips between two well-defined polarization states, where the flipping 
angle depends on a,b and on 0o, ^o- The following are two examples of switching solutions. 

As the first example we choose a = l,b = 4 with the initial values vq = e, 0o = 0, where 
e = 10^^, in which case the input laser beam is linearly polarized and the polarization state 
is close to one of the principle axes of the waveguide. Hence, the trajectory starts near the 
unstable steady states (l24l i or ( |26] |. which lie on the boundary of the red region shown in Fig.|2] 
(i). We plot V and cos ^ = cos A0 as a function of T in Fig.|6](i), showing switching behavior for 
cos -I, which is periodic and flips abruptly between the values ±1; 0, however, is an increasing 
function of T, with jumps through 2n at periodic intervals. The polarization vector experiences 
an angular flipping associated with the abrupt flipping of cosA0, however, since vq = e and 
00 = 0, the flipping angle is very small, as depicted in the inset of Fig. |6](i). Regarded as the 
trajectory of a particle of mass M in the potential V in Eq. (|32] | this motion corresponds to a 
particle moving slowly over the peaks of the potential, which are the unstable steady states, 
then sliding quickly down the valleys through the minimum values of V and back to the peaks. 
For a = 1 the potential is flat at its maximum values, since in this case V' = = V" = V'", 
hence v, are each close to zero except when moves to an adjoining maximum of V . In terms 
of the contour plots shown in Fig. Qi) this trajectory corresponds to the contour which begins 
just above the unstable steady state (orange dot) and closely follows the separatrix shown in red 
(which is the soliton solution discussed in the Appendix) with a maximum value ~ 0.4 for v. 
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Fig. 6. Switching solutions v and cos ^ = cosA(/i as functions of T for: (i) a = \,b = A 
and I'o = e, 00 = 0; (ii) a = b = 2 and vq = j,6o = e where e = 10~^. The insets show the 
polarization vectors associated with the values cos A(j) = ± 1 . 



As a second example of switching behavior we choose a~b~2 with vo = 1 /2 , 0o = £, where 
e = 10~^, which corresponds to a linearly polarized input laser beam in which the polarization 
vector makes an angle of 45° to either of the principle axes of the waveguide. Again, the initial 
value lies close to an unstable steady state jTM and a,^ lie within the red region of instability in 
Fig.|2] We plot V and cos(0/2) as functions of T in Fig.|6](ii), showing the periodicity of these 
functions and the switching behavior of cos(0/2). Since vo = 1/2, the angular flipping of the 
polarization vector is n/2, because cos(0/2) flips between values ±1, as shown in the inset of 



Fig.|6](ii). Unlike the previous example, 9 is also periodic in T with a trajectory that corresponds 
to the motion of a particle in the potential V, starting slowly near the unstable steady state (l24l 
but sliding rapidly through the potential minimum to approach an adjoining unstable steady 
state. This motion is similar to the periodic oscillations of a nonlinear pendulum (since a = b, 
see the definition of V in Eq. ( |32] |) with a large amplitude of almost 2n, and v attains nearly 
all values between 0, 1. In terms of the phase space contours shown in Fig. EJii), the motion 
corresponds to a periodic trajectory which begins near the red dot (unstable steady state) and 
again closely follows the separatrix which marks the soliton trajectory. 

5. Discussion and Conclusion 

Switching states, as defined and demonstrated here through simulation by means of a full vec- 
torial model, are attractive for practical applications, since they allow nonlinear self-flipping of 
the polarization states of light propagating in an optical waveguide. This flipping is due to the 
nonlinear interactions of the two polarizations, and has properties that depend on the total op- 
tical power and on the specific fiber parameters. These properties can in principle be employed 
to construct devices such as optical logic gates ISSll . fast optical switches and optical limiters 
Il55ll56l . in which small controlled changes in the input parameters lead to sudden changes in 
the polarization states. 

The minimum power necessary to generate such switching states is determined for any 
waveguide by the inequalities (l34l ) and, for chalcogenide optical nanowires with elliptical core 
cross sections, is summarized in Fig. |4](ii). The minimum power required in such nanowires 
is in the range 1 — lOkW which, although not practicable for CW lasers, can be achieved in 
pulsed lasers. Although we have limited our analysis to the static case, ignoring the temporal 
variation of laser light, it is still applicable to slow pulses with pulse widths in the order of 
nanoseconds depending on the dispersion of the waveguide. A more practical minimum power 
requirement that achieves switching behavior is by means of asymmetric waveguides, such as 
rib waveguides, for which Aj8 can be reduced to very small values while still having different 
field distributions for the two polarizations, as discussed in Il55ll56l . 

The nonlinear interactions of the two polarizations can be impacted by two factors that have 
not yet been investigated: (1) interactions with higher order modes in few-moded waveguides 
and, (2) contributions from nonlinear terms containing different forms of Ci . 62, i.e., nonzero 

values for the coefficients 7|iv, 7|ivi7v in Eqs. (|6][8]l. (This applies only when ei .62 is no 
longer approximately zero, as assumed in this paper). In few-moded waveguides, higher order 
modes contribute to the nonlinear phase of each polarization of the fundamental mode through 
cross phase and coherent mixing terms. Inspection of Eq. (I2) reveals that nonzero 7^^ , 7^v , 7^v 
coefficients significantly change the dynamics of nonlinear interactions of the two polarizations 
and most likely lead to different parameter regimes for the existence of periodic and solitonic 
solutions. These factors will be the subject of further studies. 

In summary, we have developed the theory of nonlinear interactions of the two polarizations 
using a full vectorial model of pulse propagation in high index subwavelength waveguides. 
This theory indicates that there is an anisotropy in the nonlinear interactions of the two polar- 
izations that originates solely from the waveguide structure. We have found all static solutions 
of the nonlinear system of equations by finding exact constants of integration, which leads to 
expressions for the general solution in terms of elliptic functions. We have analyzed the stabil- 
ity of the steady state solutions by means of a Lagrangian formalism, and have shown that there 
exist periodic switching solutions, related to a class of unstable steady states, for which there 
is an abrupt flipping of the polarization states through an angle determined by the structural 
parameters of the waveguide and the parameters of the input laser. By means of a Hamiltonian 
formalism we have analyzed all solutions, including solitons which we have shown are close to 



the switching solutions of interest. 
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Appendix 

We include here a discussion of the topological solitons which appear as solutions of Eq. ( |29] l, 
as configurations 0(t) which interpolate between the adjacent maximums of the periodic po- 
tential Y defined in Eq. ( [32] i. They define trajectories which move between adjacent unstable 
steady states with abrupt transitions, to form "kinks" which are stable against time-dependent 
perturbations. Such trajectories are visible in Fig. 0(i) and (ii) (the contours marked in red) 
as they form the separatrix between periodic solutions v, Q and nonperiodic solutions. The fact 
that solitons can occur in this way has been previously noted, see for example Chapter 9 in 
ll54l . In Fig.|5](i) the soliton is the trajectory which connects the adjacent unstable steady states 
(orange) at v = and Q ~ 0,27r,4;r. . . and similarly in Fig.|5](ii) the solitons connect the (red) 
unstable steady states. Such solutions exist on the full real line —<x> < t < oo, with appropriate 
boundary conditions, but are also solutions on any finite subset of the real line, corresponding 
to an optical fiber of finite length, with boundary values obtained from the exact solution. 

Solitons are significant in the context of switching solutions since switching behavior occurs 
precisely when solutions lie near soliton trajectories; the switching solutions shown in Fig. |6] 
(i) and (ii), for example, correspond to contours in Fig. |5](i) and (ii) which lie very close to 
the separatrix. The soliton itself is not periodic but nearby trajectories are periodic for both v 
and cos as functions of T. The abrupt transitions which characterize switching, as shown for 
example in Fig.|6] can equally be viewed as the "kinks" of a soliton, in which cos(0/2) changes 
between two distinct values over a very short T-interval, and in doing so interpolates between 
unstable steady states. We are interested here mainly in transitions between the unstable steady 
states (l24l i. since these correspond to polarization flipping, i.e. cos A0 = cos(0/2) flips between 
values ±1. There exist, however, solitons corresponding to the other unstable steady states such 
as (I26I27I I. which we also discuss briefly. 

In order to find explicit solutions, we define a potential U according to C/(0) = Vq — ^(6), 
where the shift Vo is selected such that the minimum value of U is zero. If 1 < a <2b—l, for 
example, in which case the unstable steady states (|24] | exist, we have 

Vo^l-b-^-^. (41) 

0—1 



We also define the positive "action" functional S by 



8(9, 9) 



^M{9) 9^ + U{9) 



dx. (42) 



Eqs. ( |29] l and dSTl ) follow by using Hamilton's principle of least action applied to S. We can 
write 

"V "m 



s^l \m 




(43) 



The last term takes values only on the boundary and so does not vary as 0, are varied, hence 
a local minimum of S occurs when 

(44) 




which implies MQ^ ~ 2U. Solutions of this equation, which is equivalent to Eq. (l33T l with 
c = Vq, satisfy Eqs. ( |29] | and dSTT i with the property that S < °°. Hence, for such solutions we 
have — > and 6 approaches a zero of U as |t| ^ °°. We therefore integrate Eq. (l44t or 
equivalently Eq. ( |33] ) with c — Vq. 

For the first example we select a, bin the red region in Fig. |2]for which 1 < a <2b—l, with 
c = Vo given by WH . then the soliton interpolates between the unstable steady states ( |24] |. By 
direct integration of Eq. (l33T l or (l44l we obtain 

2k- 

cos0 = l + ^^ , (45) 

1 - (K'+l)cosh v'c(t-To) 

where To is the constant of integration, and 

(fl-l)(-fl + 2/?-l) 



K 



2{b~l) 



The solution satisfies limi^i^^ cos = 1 and at T = Tq, which may be regarded as the location of 
the soliton, we have cos = — 1 . By suitable choice of sign for 0, and by choice of the branch of 
the inverse cosine function, we obtain as a function of T which either increases or decreases 
between any two adjacent zeros of the potential U at cos0 — 1. From Eq. ( |28] l we obtain v: 



1 



-b±{b- l)V'F+Tcoshv^(T-To)' 



(46) 



where the sign corresponds to either increasing or decreasing 0, and we have lim|^|^o„v(T) = 

g-l 

2(6-1)- 

As a specific example, for a ^ b = 2 and K ~ 1/2, the separatrix trajectory shown in Fig. 
|5](ii) is the parametric plot of v, as functions of the parameter t; v evidently varies between 
maximum and minimum values which occur at T = To, as can be determined directly from 
( l46l l. We can also find the solutions ( 145146b directly by solving Eq. dSTJ i. It is necessary only 
to determine Hq = //(vo,0o) by choosing vo,0o at |t| = °o, which then determines Q from 
(l38T l. For the states ( |24] | we obtain Hq = — |w, j-. and Q{v) has a repeated root at v = nfrrn'^ 
the expression ( |46] l for v may then be obtained by using the general integration formulas in 
Sections 2.266, 2.269 of Ref. l57ll . 

Solitons also exist corresponding to the unstable steady states (l25T l, provided 2b+l < a < 
— 1 and b < —I, and may be obtained from the formulas (I45]46b by means of the symmetry 
T — >■ — T, — )■ + 7r,fl — ?► —a,b -^ —b which leaves Eqs. ( 122123b invariant. The parameter K, 
for example, is now defined by k = {a + l){a — 2b— i)/2/{b + 1) which is positive in the 
orange region of Fig. |2](ii). 

Consider next the unstable states ( |26] |. which are defined only in the strip |fl| ^ 1 of the a,b 
plane. Soliton solutions take the values cos = a, v = as | t| ^ oo, and hence the Hamiltonian 
function H{v,9) defined in Eq. (l35b takes the constant value Hq — 0, which corresponds to 
c = Vo = 2(fl — b) in Eq. (1331 ). By solving v' = Q{v) we find: 

v(t) = " , (47) 

I -ab+\b — a\ cosh[vl - fl^ (t - To)] 

which exists for all |a| < 1 and b ^ a. We have limi^i^^ v(t) = and v attains its maximum 
value Vmax at T = To, with either Vfnax = {a + l)/{b+l) foib > a or else v^ax = ('5'— l)/(^— 1) 
for b < a. Having found v, we obtain cos from Eq. ( |36] | with //q = using cos 9 — {a — 



^i')-irT^^jr^T7^-r^^ -^^ = i-TT77-;rw- (49) 



bv)/{l — v), specifically 

l-fl2 

cos0(T)=fl ^ , (48) 

-fl + 77 cosh[v 1 - fl2 (t - To)] 

where rj — {b — a)/\b — a\is the sign ofb — a. We have 9 ~ a — cos0 andcos0(To) ^ —Tj. For 
the special case b ^ a with \a\ < 1, or if a = 1, we solve v^ = Q{v) directly; in the latter case 

we obtain 

2 „ , 2 

b + l + {b-l){T~To)^' 1 + (T-To)2 

As a specific example we choose a = l,b = 4, for which contour plots for constant H are 
shown in Fig. |5] (i); the (red) separatrix trajectory in particular is visible as the curve which 
connects the unstable steady states at v = 0, =0,2n This separatrix is precisely the para- 
metric plot of V, 9 given by ( |49] l. where v evidently varies between zero and its maximum value 
of 2/(i>+ 1) = 0.4 which occurs at T = Tq, while cos 9 varies between the values 1 as |t| -> °°, 
when V = 0, and — 1 at T = Tq. 

There are also solitons corresponding to the unstable steady states dZTJ i. Precise formulas can 
be obtained from Eqs. (147148b by means of the transformations 0^— 0,\'^-l — v, «—?►—« + 2fe 
which are discrete symmetries of the defining equations (I22]23b . 



